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Abstract 


The  purpose  of  this  work  is  to  compare  the  capa- 
bility of  two  different  numerical  techniques  for  pre- 
dicting the  calm-water  resistance  of  a high-speed 
sealift  trimaran.  The  vessel  is  representative  of  a 
future  advanced  naval  design  that  must  be  able  to 
travel  long  distances  with  heavy  payloads,  but  is 
constrained  geometrically  in  order  to  operate  in  aus- 
tere ports  upon  arrival.  It  has  a relatively  shallow 
draft  and  has  an  extendable  centerhull  that  pos- 
sesses a wave-piercing  bow.  This  particular  design 
has  hulls  that  are  of  essentially  equal  dimensions. 

For  comparison  of  the  codes,  a series  of  large- 
scale  model  tests  was  conducted  for  six  different 
configurations  of  the  three  hulls  (variations  of  over- 
all vessel  beam  and  length).  A thin-ship  theory  is 
shown  to  be  effective  in  predicting  resistance.  Re- 
sults for  both  the  absolute  value  and  the  relative 
change  in  resistance  due  to  change  in  the  hull  con- 
figuration, together  with  the  computational  efficacy 
of  the  method  demonstrate  utility  for  the  designer. 
The  complex  effect  on  the  resistance  due  to  the 
proximity  of  the  three  hulls  is  addressed  in  a first- 
principle  manner  to  enhance  the  prediction  from  the 
thin-ship  theory. 

The  commercial  computational  fluid  dynamics 
code  FLUENT  is  shown  to  more  accurately  predict 
the  resistance  when  the  calculations  are  performed 
with  the  body  fixed  at  the  experimental  sinkage  and 
trim.  The  temporal  expense  associated  with  com- 
putational fluid  dynamics  is  reduced  in  this  study 
by  using  very  simple  computational  gilds.  Also,  the 
solution  convergence  is  enhanced  by  using  polyhe- 
dral finite  volumes  which  improve  cell  skewness. 


1 Introduction 


1.1  Background 

The  trimaran  vessel  has  been  under  consideration 
for  many  years  now.  However,  it  is  only  in  re- 
cent years  that  the  concept  has  been  subjected  to 
high-level  technical  scrutiny.  The  reader  is  referred 
to  books  by  Dubrovsky  and  Lyakhovitsky  (2001) 
and  Dubrovsky  (2004),  where  the  general  advan- 
tages and  disadvantages  are  explained.  Other  pub- 
lications of  a design  or  descriptive  nature  include 
those  of  Summers  and  Eddison  (1995),  Pattison 
and  Zhang  (1995),  van  Griethuysen,  Bucknall,  and 
Zhang  (2001).  In  these  papers,  the  principal  inter- 
est, has  been  in  the  application  of  the  concept  to 
naval  vessels,  so  that  the  emphasis  has  been  placed 
on  matters  of  layout  and  operability,  rather  than  de- 
tails of  resistance  and  structural  strength.  The  pre- 
cise categorization  of  vessel  type  is  often  difficult; 
an  interesting  example  of  a hybrid  form,  namely  a 
small-waterplane-area  trimaran,  was  described  by 
Lamb  (2004).  A recent  example  of  a successful  civil 
application  of  the  trimaran  is  that  detailed  by  Arm- 
strong (2004).  In  the  development  of  that  vessel, 
effort  was  invested  in  determining  the  optimal  geo- 
metric configuration  to  minimize  both  overall  resis- 
tance and  motions  in  a seaway. 

A fundamental  question,  which  must  be  asked, 
relates  to  the  potential  resistance  advantage  of  a 
trimaran  over  more  traditional  hullforms,  such  as 
monohulls  and  catamarans.  Over  three  decades  ago, 
Narita  (1976)  examined  this  problem,  using  a multi- 
body extension  to  the  hallmark  linearized  theory  of 
Michell  (1898)  for  a monohull.  This  matter  was  also 


studied  in  depth  by  Doctors  (1999),  who  compared 
both  the  wave  resistance  and  the  total  resistance  of 
vessels  possessing  between  one  and  six  subhulls. 

One  is  strongly  interested  in  minimizing  wave 
resistance  if  the  concern  is  about  the  wave  genera- 
tion, such  as  for  river-based  ferries,  when  bank  ero- 
sion must  be  avoided.  In  this  case,  increasing  the 
number  of  subhulls  is  generally  favorable,  although 
there  is  little  additional  reduction  in  wave  genera- 
tion after  two  or  three  subhulls  — depending  on  the 
speed  range  of  interest.  If  there  is  a primary  con- 
cern to  minimize  total  resistance,  it  is  found  that 
the  increase  in  wetted  surface  area  is  a major  draw- 
back of  a vessel  with  many  subhulls.  Indeed,  if  a 
trimaran  is  to  compete  on  this  basis,  it  is  vital  to 
keep  to  a minimum  the  size  of  the  sidehulls.  In  such 
circumstances,  the  sidehulls  may  be  referred  to  as 
outriggers. 

In  this  current  research,  we  also  principally 
wish  to  investigate  the  resistance  characteristics  of 
trimarans,  because  the  application  is  to  sealift  ves- 
sels, in  which  the  operational  range  is  to  be  max- 
imized. Thus,  there  is  a need  to  minimize  the  to- 
tal resistance.  Clearly,  the  most  obvious  idea  is  to 
consider  different  longitudinal  positions  (the  stag- 
ger) and  different  lateral  positions  (the  spacing) 
of  the  sidehulls.  In  this  regard,  one  can  refer  to 
the  work  of  Wood  (1988),  Wilson  and  Hsu  (1992), 
Gale,  Hall  and  Hartley  (1996),  Ackers,  Michael,  Tre- 
dennick,  Landen,  Miller  III,  Sodowsky  and  Hadler 
(1997),  Battistin,  Danielli  and  Zotti  (2000),  Doc- 
tors and  Scrace  (2003),  Day,  Clelland  and  Nixon 
(2003),  Begovic,  Bove,  Bruzzone,  Caldarella,  Cas- 
sella,  Ferrando,  Tincani  and  Zotti  (2005),  and  Ya- 
sukawa  (2005).  A common  theme  is  that  the  side- 
hulls  should  be  located  aft  in  order  to  reduce  the 
resistance.  Obviously,  the  precise  answer  must  de- 
pend on  the  relative  displacements  of  the  centerhull 
and  the  sidehull,  as  well  as  the  intended  operational 
Froude  number. 

The  specific  matter  of  optimization  of  the  con- 
figuration of  the  trimaran  is  closely  linked  to  the 
previous  comments.  This  topic  has  been  investi- 
gated by  a number  of  researchers,  including  Suzuki 
and  Ikehata  (1993),  Yang,  Soto,  Lohner  and  No- 
blesse (2002),  Yang,  Noblesse  and  Lohner  (2002), 
Degiuli,  Werner  and  Doliner  (2003),  and  Brizzolara, 
Bruzzone,  and  Tincani  (2005).  When  studied  in  a 
theoretical  or  computational  sense,  it  must  be  un- 
derstood that  the  problem  is  nonlinear.  That  is,  the 
only  reliable  approach  to  determining  the  global  op- 


timum is  through  a search  technique  that  is  coded 
in  order  to  avoid  local  false  optimums. 

It  is  necessary  to  test  theoretical  prediction 
methods  against  experiments.  This  is  most  conve- 
niently achieved  at  model  scale  in  a towing  tank. 
Examples  of  such  effort  is  characterized  by  the 
papers  by  Scrace  (2000),  Kennell  (2004),  Mizine, 
Amromin,  Crook,  Day  and  Korpus  (2004),  Colic- 
ehio,  Colagrossi,  Lugni  and  Faltinsen  (2005),  and 
Degiuli,  Werner  and  Zotti  (2005).  In  general,  it 
can  be  stated  that  current  methods  based  on  the 
Michell  (1898)  theory  or  computational  fluid  dy- 
namics (CFD)  provide  useful  results.  Nevertheless, 
the  predictions  are  often  not  as  accurate  as  those 
for  catamarans.  It  is  believed  that  this  is  due  to 
the  possible  strong  flow  interferences  between  the 
sidehulls  and  the  centerhull,  which  can  create  an 
adverse  effect  (that  is,  a large  increase)  on  the  fric- 
tional form  factor.  We  will  specifically  refer  to  this 
point  again  later  in  this  paper. 


1.2  Current  Work 

In  the  current  work,  we  describe  an  exhaustive  in- 
vestigation into  the  resistance  for  a candidate  high- 
speed sealift  trimaran.  The  purpose  of  the  effort 
was  twofold:  to  consider  a suitable  trimaran  design 
and  to  test,  our  available  computational  tools. 

The  software  developed  by  Doctors  (1999)  and 
Doctors  and  Scrace  (2003),  based  on  extensions  to 
the  Michell  (1898)  theory,  represented  one  focus  of 
the  computational  work.  Computations  were  per- 
formed in  a “blind  fashion”,  in  that  the  experimen- 
tal results  for  the  resistance  were  not  initially  made 
available  to  those  persons  performing  the  calcula- 
tions. Of  course,  the  opportunity  was  later  provided 
to  learn  from  the  comparisons  of  theory  and  exper- 
iment and,  consequently,  to  improve  the  software. 

The  second  focus  of  the  work  was  based  on  a 
Reynolds-Averaged  Navier-Stokes  (RANS)  compu- 
tational fluid  dynamics  (CFD)  approach  using  the 
commercially  available  software  FLUENT.  The  ini- 
tial intentions  were  to  perform  the  computations  in 
a “blind  fashion”,  as  in  the  Doctors  results;  how- 
ever, difficulties  were  encountered  in  dynamically 
predicting  the  sinkage  and  trim  of  the  model  as  part 
of  the  computations.  Therefore,  several  fixed  model 
conditions  were  simulated,  using  the  experimentally 
measured  sinkage  and  trim  for  the  model.  It  is  not 


Table  1,  for  the  loading  cases  where  the  draft  is  set 
at  6.510  m. 


(a)  Pictorial  of  Configuration  6 


• / i /////} 

\:/m 


(b)  Side  Hull 


(c)  Center  Hull 


Figure  1:  Vessel  Geometry 


yet  clear  what  the  limiting  factor  is  for  predicting 
the  sinkage  and  trim  of  the  model.  Previous  in- 
vestigations by  some  of  the  authors  have  success- 
fully computed  ship  model  behaviors  for  pitch  and 
heave  in  the  presence  of  incident  waves  (Sasanapuri, 
Shirodkar,  Wilson,  Kadam,  and  Rhee,  2007),  and 
further  studies  are  planned  to  include  the  dynamic 
motion  of  the  model  as  part  of  the  computations. 


2.2  Advantages  of  Concept 

It  is  well  known  that  a resistance-optimal  conven- 
tional trimaran  is  one  in  which  the  sidehulls  possess 
a very  small  displacement  compared  to  that  of  the 
centerhull.  However,  in  the  present  application,  it 
is  intended  that  the  trimaran  will  be  operated  in 
essentially  two  different  (and  diverse)  modes. 

In  the  cruising  mode,  the  sidehulls  are  to  be 
staggered  rearward  by  up  to  50%  of  the  length  of 
the  subhulls.  In  this  configuration,  the  great  over- 
all length  of  the  vessel  renders  it  somewhat  simi- 
lar to  a slender  monohull  from  a resistance  point  of 
view.  Thus,  its  resistance  qualities  are  expected  to 
be  good.  In  the  arrival  or  departure  mode,  the  cen- 
terhull  can  be  retracted  so  that  it  is  aligned  with  the 
sidehulls.  While  the  resistance  is  certain  to  be  un- 
reasonably high  in  this  configuration,  it  will  permit 
the  vessel  to  be  docked  in  an  austere  port. 

This  vessel  was  designed  by  Dr  R.  Scher,  of 
Alion  Science  and  Technology,  in  Alexandria,  Vir- 
ginia. The  designing  of  this  vessel  formed  part  of  the 
project  Architectural  Concepts  and  Hydrodynamic 
Technologies  for  High  Speed  Sealift  to  Austere  Ports: 
Subtopic  B:  Computational  Approach  and  Hydrody- 
namic Tools,  supported  by  the  US  Office  of  Naval 
Research,  under  the  supervision  of  Project  Manager 
Dr.  L.  Patrick  Purtell. 


2 Design  of  Candidate 

2.1  Geometry 

The  vessel  chosen  for  our  investigation  is  a trimaran, 
in  which  the  sidehulls  and  the  centerhull  are  essen- 
tially equal  in  length,  beam,  draft,  and  displace- 
ment. The  vessel  is  depicted  in  Figure  1,  where  the 
perspective  view  of  the  extended  vessel  can  be  seen 
along  with  the  body  plans  of  the  side  and  center 
hulls. 

Both  the  centerhull  and  the  sidehulls  possess 
a transom  stern,  which  would  permit  the  installa- 
tion of  suitable  waterjets.  The  principal  geometric 
and  hydrostatic  data  for  the  subhulls  are  listed  in 


3 Towing- Tank  Tests 

3.1  Physical  Description  of  Model 

The  bare-hull  model  resistance  tests  were  conducted 
at  the  Maritime  Research  Institute  Netherlands — 
MARIN.  Three  individual  hulls  were  constructed 
out  of  wood,  with  a scale  factor  of  1:34,  and  con- 
nected together  using  aluminum  beams  which  al- 
lowed for  easy  reconfiguration.  Figure  2 shows  an 
image  of  the  model  in  Configuration  4. 

The  model  was  attached  to  the  carriage  and 
towed  with  a vertical  heave  staff  and  gimbal  so  that 
the  model  was  free  to  heave  and  pitch.  The  longi- 


Table  1:  Nominal  Particulars  of  HSSL  ALT  Subhulls 


Item 

Symbol 

Units 

Center- 

hull 

Side- 

hull 

Displacement  mass 

A 

t 

5827 

5964 

Waterline  length 

L 

m 

168.1 

168.6 

Waterline  beam 

B 

m 

10.03 

10.04 

Draft 

T 

m 

6.510 

6.510 

Waterplane-area  coefficient 

CWp 

0.7770 

0.8091 

Maximum  section  coefficient 

CM 

0.7860 

0.7863 

Block  coefficient 

Cb 

0.5174 

0.5279 

Prismatic  coefficient 

cP 

0.6583 

0.6714 

Slenderness  coefficient 

L/V1/3 

9.424 

9.375 

Figure  2:  Image  of  model  in  Configuration  4 


tudinal  location  of  the  tow  point  corresponded  to 
the  longitudinal  center  of  buoyancy  of  the  trimaran 
(which  varied  depending  on  the  stagger  of  the  side- 
hulls)  . 


3.2  Test  Configurations 

The  towing-tank  model  was  operated  in  six  differ- 
ent loading  configurations,  as  listed  in  Table  2 and 
viewed  in  Figure  3.  Essentially,  this  consisted  of 
three  different  staggers  of  the  sidehulls  and  two  dif- 
ferent offsets  S2/2  of  the  eenterplane  of  the  sidehulls 
from  the  eenterplane  of  the  centerhull.  During  the 
planning  phase  of  the  experiments,  we  had  intended 
that  the  draft  should  be  identical  for  all  test  con- 
figurations. Unfortunately,  the  freeboard  proved  to 
be  too  small  for  the  reduced  dimensions  for  Config- 
uration 4 and  Configuration  5;  as  a result,  an  exces- 
sive amount  of  spray  impacted  on  the  bridging  cross 
structure  connecting  the  three  subhulls.  This  defi- 
ciency in  the  design  of  the  trimaran  could  easily  be 
corrected  in  any  planned  extension  to  the  project. 


~ " — — 

(b)  Configuration  2 


(c)  Configuration  3 


(d)  Configuration  4 


(e)  Configuration  5 


(f)  Configuration  6 


Figure  3:  Bottom  View  of  the  6 Hull  Configurations 


It  is  recommended  that  the  wetdeck  height, 
which  in  the  current  design  is  10.63  m,  be  increased 


by  a minimum  of  2 m.  This  is  based  on  observations 
of  the  resistance  tests  and  of  the  seakeeping  exper- 
iments that  were  performed  on  the  model.  (Sea- 
keeping experiments  were  also  done  solely  for  Con- 
figuration 6 and  the  model  was  self-propelled  with 
a single  water  jet  fitted  in  each  hull.)  Two  factors 
contributed  to  the  relative  rise  of  the  free  surface  on 
body.  Firstly,  the  geometry  of  the  shorter  Configu- 
rations 1,  2,  4 and  5 places  the  forward  portion  of 
the  sidehulls  in  the  peak  of  the  dominant-diverging 
wave  of  the  centerhull.  Secondly,  the  vessel  sinks 
downward  and  trims  bow  down  at  a speed  of  ap- 
proximately 30  knots.  Figures  4 and  5 depict  the 
effect  of  forward  speed  on  the  running  attitude  of 
the  body,  and  thus  the  freeboard.  Figure  4 shows 
in  the  deep  draft  Configuration  4 at  30  knots,  and 
the  vanishing  freeboard  can  be  seen.  Note  that  the 
models  were  tested  with  the  portion  of  the  geome- 
try from  the  baseline  to  the  wetdeck  plane.  Figure  5 
shows  how  the  vessel  attitude  changes  at  a higher 
speed.  The  bow-up  trim  and  rise  in  the  vessel  act 
to  increase  the  freeboard  on  the  vessel. 

Other  observations  from  the  seakeeping  exper- 
iments suggest  that  a 2 m rise  in  the  wetdeck  should 
be  sufficient  to  allow  clearance  for  the  spray  and  jet 
of  water  exiting  the  centerhull.  Finally,  a rise  in  the 
wetdeck  will  allow  for  larger  free-board  on  the  cen- 
terhull, which  is  currently  limited  by  the  fact  that  it 
must  retract  underneath  the  sidehulls  for  entry  into 
Port  Austere. 


Figure  5:  Image  of  model  underway,  in  Configura- 
tion 1,  at  a speed  of  44  knots 

4 Theoretical  Approach 


4.1  Thin-Ship  Analysis  of  a 
Transom-Stern  Multihull 


Figure  4:  Image  of  model  underway,  in  the  heavy 
Configuration  4 (ballasted  to  a draft  of  6.51  m),  at 
a,  speed  of  30  knots 


A computational  approach  based  on  the  work  of 
Doctors  and  Day  (1997)  was  used  to  predict  the  re- 
sistance of  the  1/34-scale  model.  This  approach  has 
its  roots  in  the  linearized  theory  of  Michell  (1898). 
The  current  version  of  the  software  contains  a num- 
ber of  enhancements. 

Firstly,  the  program  can  handle  any  number  of 
subhulls  up  to  six,  by  employing  additional  center- 
plane  source  distributions.  The  transverse  velocities 
induced  by  one  subhull  should  be  corrected  for  by 
the  use  of  centerplane  dipole  distributions;  however, 
this  is  ignored  for  the  sake  of  considerable  simplifi- 
cation to  the  coding  of  the  software. 

Secondly,  an  enhanced  model  for  the  transom- 
stern  flow  has  been  under  development  for  ten  years 
now.  The  current  model  allows  for  a realistic  hy- 
drodynamic addition  to  the  vessel  by  means  of  a 
virtual  extension  to  the  vessel  in  order  to  represent 
the  presence  of  the  transom-stern  hollow.  This  has 


been  described  by  Doctors  and  Beck  (2005),  Doctors 
(2006a  and  2006b),  and  Maki,  Doctors,  Beck,  and 
Troesch  (2006).  In  this  way,  the  effective  hydrody- 
namic length  of  the  subhulls  grows  with  increasing 
vessel  speed  in  a physically  realistic  manner  that  has 
been  substantiated  by  numerous  towing-tank  exper- 
iments. 

Thirdly,  the  transom-stern  model  includes  a 
procedure  for  estimating  the  progressive  unwetting 
or  ventilation  of  the  transom.  This  feature  permits 
a practical  estimate  of  the  transom-stern  or  hydro- 
static drag  suffered  by  such  vessels. 

Fourthly,  as  a special  task  inspired  by  the 
model  experiments  to  be  described  here,  a further 
procedure  was  incorporated  in  the  computer  code 
in  order  to  take  into  account  the  viscous  interac- 
tion between  the  subhulls.  This  simplistic  proce- 
dure makes  an  allowance  that  the  subhulls  are  in 
the  “shadow”  of  each  other  to  some  degree,  as  seen 
in  a profile  or  a side  view.  This  shadow  will  be  great- 
est when  the  subhulls  are  abreast  of  each  other  and 
will  be  particularly  strong  when  they  are  of  similar 
size.  In  the  present  theory,  the  local  gap  between 
the  subhulls  below  the  loaded  waterplane  was  com- 
puted. From  this,  an  estimate  of  the  increase  in  lo- 
cal water  speed  was  obtained.  This  increase  in  speed 
was  applied  to  the  usual  formula  for  frictional  resis- 
tance, thus  leading  to  an  approximate  technique  for 
determining  the  frictional  form  factor. 

Of  course,  it  is  clearly  understood  that  a num- 
ber of  physical  phenomena  is  neglected  here.  Ne- 
glected effects  include  the  deformed  shape  of  the 
free  surface  and  the  fact  that  the  flow  would  tend 
to  skirt  around  the  outside  of  the  vessel  in  cases 
where  the  intersubhull  gap  is  very  small.  However, 
we  wished  to  preserve  the  characteristically  simple 
approach  of  the  software  as  well  as  its  fast  compu- 
tational behavior. 


4.2  Analysis  Using  RANS 

Computational  fluid  dynamics  codes  have  demon- 
strated increasing  fidelity  in  recent  years  in  predict- 
ing the  Kelvin  wave  structure  for  a variety  of  ships 
and  marine  craft.  Advancements  in  computational 
methods  for  free-surfaee  predictions  using  level-set 
and  volume-of-fluid  (VOF)  techniques  have  demon- 
strated marked  improvements  in  the  ability  to  ac- 
curately predict  flows  around  surface  ships  that  are 


moving  at  speeds  that  cause  breaking  waves.  Re- 
cent advances  in  computing  technologies  and,  in 
particular,  parallel-processing  techniques,  have  also 
greatly  enhanced  the  ability  to  perform  simulations 
with  increased  spatial  resolution  of  the  flow  field 
and  with  faster  compute  times.  These  improve- 
ments, along  with  growing  user  experience  in  ap- 
plying RANS  methods  to  these  types  of  problems, 
are  beginning  to  demonstrate  the  ability  of  CFD 
methods  to  impact  design  and  evaluation  for  ma- 
rine hydrodynamic  problems.  A review  of  applying 
RANS  techniques  to  a variety  of  surface  ship  con- 
figurations was  given  by  Gorski  (2004).  An  assess- 
ment of  different  computational  tools  was  given  in 
Wilson,  Fu,  Fullerton,  and  Gorski  (2006). 

The  commercially  available  viscous  RANS  flow 
solver  FLUENT  was  used  to  predict  the  resistance 
for  the  1 :34-scale  model.  The  computations  focused 
on  Configuration  6 (Figure  3(f))  as  this  was  ex- 
pected to  be  the  best  performing  design  configu- 
ration for  resistance  and  seakeeping.  The  compu- 
tations were  performed  using  FLUENT  V6.3  and 
the  parallel  computational  resources  at  the  Ship  En- 
gineering and  Analysis  Technology  Center,  located 
at  the  Naval  Surface  Warfare  Center,  Carderock 
Division  (NSWCCD)  and  the  Aeronautical  Sys- 
tems Center  (ASC)  Major  Shared  Resource  Center 
(MSRC)  High  Performance  Computing  facility. 

The  CFD  software  uses  a cell-centered  finite- 
volume  method  and  allows  a variety  of  differ- 
ent computational  element  (cell)  types,  including 
quadrilateral,  tetrahedral,  hexahedral,  pyramidal, 
prismatic,  and  hybrid  meshes.  The  latest  version 
of  the  software  also  includes  the  ability  to  con- 
vert tetrahedral  mesh  zones  into  polyhedra,  which 
are  constructed  from  an  arbitrary  number  of  sides. 
Computational  meshes  constructed  from  polyhedra 
have  several  advantages,  including  reduced  compu- 
tational cell  count  and  greatly  improved  mesh  qual- 
ity through  reduced  cell  skewness.  Both  of  these 
attributes  have  a significant  positive  effect  on  solu- 
tion convergence. 

The  computational  grids  for  performing  the 
CFD  calculations  employed  hexahedral  volumes  in 
the  far-field,  boundary-layer  prisms  on  the  hulls, 
and  polyhedral  volumes  to  connect  the  boundary- 
layer  prism  caps  to  the  far-field  volume.  This  de- 
composition exploits  the  advantages  of  each  of  the 
different  types  of  finite- volume  cell,  while  also  con- 
sidering the  desire  to  minimize  the  difficulty  in  grid 
generation. 


Table  2:  Details  of  Model  Experiments 


Config- 

uration 

Overall 
Length 
L (m) 

Overall 
Beam 
B (m) 

Sidehull 

Stagger 
r2  (m) 

Sidehull 
Offset 
s2/2  (m) 

Draft 
T (m) 

Displace- 

ment 

A (t) 

Speed 
Range 
U (kn) 

1 

168.6 

56.0 

0.0 

23.0 

6.51 

17,756 

20-56 

2 

210.6 

56.0 

-42.5 

23.0 

6.51 

17,756 

20-60 

3 

253.1 

56.0 

-85.0 

23.0 

6.51 

17,756 

20-54 

4 

168.2 

39.9 

0.0 

15.0 

5.19 

12,341 

20-54 

5 

210.2 

39.9 

-42.5 

15.0 

5.19 

12,341 

20-48 

6 

253.1 

40.0 

-85.0 

15.0 

6.51 

17,756 

20-58 

Figure  6 shows  the  surface  mesh  on  the  center- 
hull  and  centerplane  where  each  of  the  three  differ- 
ent cell  types  can  be  seen.  Also  visible  in  this  Figure 
is  the  thickness  over  which  the  free  surface  is  arti- 
ficially spread  due  to  numerical  discretization.  The 
thickness  is  seen  to  be  at  most  three  cells,  a value 
that  it  not  intolerable  because  even  the  most  ad- 
vanced and  thus  complex  advection  schemes  rarely 
compress  the  interface  to  fewer  than  two  cells. 

The  hexahedral  volumes  are  effective  in  solv- 
ing both  the  undisturbed  free-surface  upstream  of 
the  body  and  the  diffracted  waves  away  from  the 
body,  with  a minimum  number  of  finite  volumes. 
This  is  contrary  to  the  polyhedra  which  are  better 
suited  for  filling  volumes  of  the  flow  domain  that 
are  defined  by  complex  geometry  but  compromise 
interface  resolution. 

The  polyhedra  are  touted  to  reduce  the  cell 
count  by  up  to  a factor  of  five  when  compared  to 
a tetrahedral  unstructured  grid,  but  such  a large 
reduction  of  unknowns  must  have  an  impact  on  the 
resolution  of  the  flow  field.  In  this  application,  the 
polyhedra  were  created  from  a tetrahedral  mesh, 
and  as  to  not  excessively  reduce  the  resolution  by 
losing  too  many  cells,  the  original  tetrahedral  mesh 
was  created  with  a much  greater  density. 

The  division  of  two  separate  domains  per- 
mits that  one  far-held  grid  be  generated  and  used 
for  all  calculations  and  a new  inner  domain  be 
generated  for  each  simulation  that  requires  the 
hull  be  arranged  with  a different  sinkage  and 
trim.  The  outer  domain  had  extents  of  ( x , ?/,  z)  £ 
(-20.0,  0.0,  -20.5)  : (40.0, 20.0,  20.5)  in  model  scale 
units  of  meter.  The  inner  domain  extended  from 
(x,y,z)  £ (-8.0,  0.0, -0.5)  : (1.0,  .75, 0.5).  The 
Cartesian  coordinate  system  is  situated  at  the  cen- 
terline, transom,  calmwater  intersection,  and  ori- 


ented with  z vertically  upwards  and  x downstream. 
A layer  of  boundary-layer  prisms  was  generated  on 
each  hull.  The  first  cell  had  a height  of  0.001  m and 
at  the  speed  of  40  kn  this  resulted  in  wall  y+  values 
between  2.0  and  100,  with  the  average  value  of  39. 
The  grids  used  in  this  work  had  a total  cell  count  of 
approximately  1.6M,  composed  of  1.34M  polyhedra 
and  265K  hexahedra. 

In  the  present  study,  the  convective  terms 
were  discretized  using  the  QUICK  scheme  (Leonard, 
1979).  This  scheme  uses  a weighted  average  of  up- 
wind and  central  second-order  interpolation  meth- 
ods, which  is  typically  more  accurate  for  structured 
grids  aligned  with  the  flow  direction.  For  unstruc- 
tured meshes,  the  scheme  reverts  to  the  second- 
order  upwind  method.  Therefore,  this  scheme  is 
very  appropriate  for  the  present  case,  which  involves 
a hybrid  mesh  that  contains  both  structured  hexa- 
hedral elements  which  are  mostly  aligned  with  the 
general  flow  direction,  as  well  as  polyhedra  elements 
for  flexibility  in  resolving  the  model  geometry  and 
near-body  flow.  The  discretization  of  the  volume 
fraction  equation  used  a modified  version  of  the 
high-resolution  interface  capturing  (HRIC)  method 
(Muzaferija,  Peric,  Sames,  and  Schellin,  1998).  The 
HRIC  scheme  has  been  shown  to  be  more  accurate 
than  the  QUICK  or  other  second-order  discretiza- 
tion schemes,  and  is  less  computationally  expensive 
than  using  a complete  geometric  reconstruction  of 
the  interface.  The  turbulence  closure  was  accom- 
plished using  a model  that  is  a variant  of  the  k — u> 
model  as  described  in  Wilcox  (1998). 

The  simulations  were  performed  as  fully  un- 
steady in  time,  marching  to  an  equilibrium  solu- 
tion. The  time  advancement  was  accomplished  us- 
ing second-order  backward  differencing,  and  alge- 
braic multi-grid  methods  were  used  to  aid  in  the 
solution  convergence. 


Table  3:  Nomenclature 


Figure  6:  CFD  grid  and  flow  solution.  Contours 
of  volume  fraction  depict  the  location  of  the  free- 
surface  for  the  40  knot  speed  case 


Symbol 

Meaning 

U 

Speed  of  vessel 

W 

Displacement  weight 

fs 

Surface- velocity  limit  factor 

r2 

Stagger  of  sidehull 

s2/2 

Offset  of  sidehull 

A 

Displacement  mass 

Rp 

Frictional  resistance 

Rh 

Hydrostatic  resistance 

Rt 

Total  resistance 

Rw 

Wave  resistance 

^air 

Air  resistance 

ted  with  the  velocity  expanded  to  full  scale  using 
Froude  scaling. 

We  first  consider  Configuration  1 in  Fig- 
ure 7(a).  The  wave  resistance  Rw  possesses  its  typ- 
ical maximum  at  the  “hump”  condition.  The  hydro- 
static resistance  Rh  is  seen  to  be  relatively  small, 
indicating  that  little  of  the  transom  is  immersed  in 
the  water.  The  frictional  resistance  Rp  is  estimated 
from  the  1957 ITTC  line  and  is  seen  to  be  the  major 
drag  component  at  the  high-speed  end  of  the  speed 
range,  namely  at  50  to  60  knots.  The  simple  sum 
of  these  components  is  indicated  by  the  first  of  the 
curves  for  total  resistance  Rp.  This  predicted  re- 
sult, may  be  compared  with  the  experimental  data 
as  represented  by  the  symbols. 


5 Results 

5.1  Comparison  of  Theory  and  Ex- 
periment 


It  can  be  seen  that  this  first  prediction  of  the 
total  resistance,  indicated  by  the  symbol  fa  = 1, 
falls  short  of  the  experimental  resistance.  This  dis- 
crepancy increases  with  speed.  It  is  believed  that 
most  of  this  discrepancy  can  be  traced  back  to 
the  interactions  between  the  subhulls  increasing  the 
frictional  drag,  as  noted  earlier  in  this  paper. 


We  turn  to  Figure  7,  whose  six  parts  show  the  ex- 
perimental and  computed  results  for  the  above  men- 
tioned six  configurations.  The  symbols  on  these 
plots  are  listed  in  Table  3. 

The  curves  for  the  resistance  components  are 
all  represented  in  a dimensionless  form  known  as 
the  specific  resistance,  namely,  the  ratio  of  the  re- 
sistance to  the  vessel  weight.  It  is  important  to  clar- 
ify that  the  calculations  are  all  executed  at  model 
scale  in  order  to  be  directly  comparable  with  the 
model  experiments.  However,  the  results  are  plot- 


To  emphasize  this  point,  one  has  only  to  note 
that  the  subhulls  each  possess  a maximum  local 
beam  of  approximately  10  m at  the  waterline.  With 
an  overall  sidehull-centerplane  spacing  S2  of  46  m, 
this  leads  to  a minimum  gap  of  13  m.  This  suggests 
that  the  water  has  been  channeled  from  an  initial 
width  of  23  m down  to  13  m,  implying  a large  in- 
crease in  the  local  speed  of  the  water  over  the  sur- 
face of  the  hull.  On  the  other  hand,  it  should  be 
noted  that  this  effect  will  be  much  less  on  other 
parts  of  the  subhull  surface.  The  current  estimate 
of  the  velocity  increase  does  take  this  effect  fully 


Figure  7:  Components  of  Resistance  at  Model  Scale,  (o)  experiment;  (— ) thin-ship  theory;  (□)  FLUENT 
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Figure  8:  Resistance  coefficients  at  model  scale, 
Configuration  6.  (o)  experiment;  (— ) thin-ship  the- 
ory with  no  form  factor;  (□)  FLUENT 


into  account  in  the  third  theoretical  prediction,  in- 
dicated by  fs  = oo. 

It  is  considered  that  the  prescription  explained 
here  overestimates  the  effect  under  study,  as  noted 
when  the  vessel  speed  U exceeds  44  knots.  Thus, 
yet  another  curve  is  shown,  indicated  by  fs  = 1.2. 
The  meaning  of  this  is  that  the  coding  places  an 
upper  limit  on  the  subhull  surface  velocity  of  1.2 
times  the  vessel  speed. 

The  other  five  parts  of  Figure  7 pertain  to  the 
other  five  configurations.  On  the  whole,  it  is  con- 
firmed that  the  surface- velocity  limit  factor  fs  = 1.2 
does  seem  to  be  a reasonable  compromise  for  achiev- 
ing the  best  correlation  between  theory  and  exper- 
iment. 

Figure  7(f)  also  shows  some  calculations  based 
on  the  FLUENT  software.  Because  computational 
times  are  much  greater  for  these  very  extensive  and 
detailed  CFD  studies,  a limited  number  of  com- 
puted points  only  is  depicted  here.  It  is  indeed  most 
encouraging  to  observe  the  excellent  matching  be- 
tween the  FLUENT  data  points  and  the  experimen- 
tal data. 

Figure  8 contains  the  data  of  Figure  7(f)  re- 
plotted non-dimensionally  as  the  resistance  coeffi- 
cient, Ct  = Rt/{^/2pU2S),  at  model  scale.  Thus, 
there  is  no  scaling  included  in  the  experimental  re- 
sult. The  thin-ship  theory  results  are  calculated 
with  the  surface- velocity  limit  factor  fs  = 1.0.  That 


is,  a frictional  form  factor  of  1.0  is  employed  here. 
Clearly,  the  use  a larger  value  of  the  form  factor, 
such  as  1 + k = 1.20,  would  provide  excellent  pre- 
dictions for  the  thin-ship  theory  code.  The  results 
show  good  agreement  by  both  codes  with  the  experi- 
mental data,  but  more  accurate  results  are  provided 
by  the  FLUENT  RANS  code.  The  average  relative 
difference  between  prediction  and  experiment,  when 
using  the  speed  range  of  20  to  50  knots,  is  9.9%  and 
1.8%  for  the  thin-ship  theory  and  CFD  results  re- 
spectively. When  comparing  the  two  computational 
methods,  we  must  recall  the  extra  expense  of  the  in- 
creased accuracy.  The  time  required  to  complete  a 
converged  drag  result  for  a single  forward  speed  is 
approximately  24  to  48  hours  using  24  processors 
for  the  RANS  method  and  only  about  1 second  on 
a single  processor  for  the  thin-ship  theory  code. 

A sample  prediction  of  free-surface  elevation 
from  the  RANS  computations  for  an  equivalent  for- 
ward speed  of  40  knots  is  shown  in  Figures  9.  In- 
terrogation of  the  three-dimensional  flow  solution  is 
very  useful  in  understanding  the  resulting  Kelvin 
wake  pattern,  as  well  as  the  stern  wake  flows  and 
the  interference  effects  between  the  center  and  side 
hulls.  A closer  examination  of  these  figures  also 
points  out  how  the  solution  of  the  wave  profiles  be- 
comes somewhat  degraded  as  the  flow  information 
is  transmitted  across  the  non-conformal  grid  inter- 
face separating  the  near-body  region  from  the  outer 
domain.  Because  the  focus  of  this  effort  was  to  ex- 
amine the  ability  of  the  tools  to  predict  the  model 
resistance,  this  was  deemed  an  acceptable  way  to 
relax  the  mesh  resolution  requirements,  while  not 
significantly  impacting  the  resistance  computation, 
as  evidenced  by  the  results  in  Figures  7(f)  and  8.  Of 
course,  if  one  were  interested  in  accurately  predict- 
ing the  wave  behavior  further  away  from  the  hull, 
then  commensurate  mesh  resolution  would  be  re- 
quired. Wilson,  Fu,  Fullerton,  and  Gorski  (2006) 
have  shown  that  increased  grid  resolution  enables  a 
more  accurate  prediction  of  the  wave  elevation,  but 
the  effect  on  resistance  prediction  is  still  unclear. 

It  is  also  conjectured  that  there  is  perhaps  a 
fundamental  physical  effect  of  the  sidehull  wave  in- 
terference drag  which  is  being  better  represented  by 
the  RANS  method.  The  sidehull  wave  interference 
effect  can  be  seen  in  Figure  10.  As  pointed  out 
earlier,  this  effect  increases  with  increasing  forward 
speed;  hence,  the  discrepancy  between  the  two  com- 
putational methods  grows  larger  at  the  higher  speed 
range.  We  restate  that  the  RANS  simulations  are 
conducted  without  including  the  sinkage  and  trim 


Figure  9:  Free-surface  contours  colored  by  elevation 


prediction  as  part  of  the  computation,  but  rather 
the  model  attitude  is  fixed  based  on  the  experimen- 
tally measured  sinkage  and  trim.  This  most  likely 
provides  a reduction  in  the  error  by  eliminating  the 
difference  between  the  predicted  and  physical  model 
sinkage  and  trim  that  would  propagate  into  errors 
in  the  predicted  hull  forces. 


5.2  Effect  of  Changes  to  Configura- 
tion 


We  now  proceed  to  Figure  11,  in  which  we  compare 
the  change  in  total  resistance  when  comparing  two 
different  configurations. 

Figure  11(a)  shows  the  fractional  change  in 
specific  total  resistance  A Rt/W  between  Config- 
uration 1 and  either  Configuration  2 or  Configura- 
tion 3.  We  can  point  out  two  important  features. 
Firstly,  the  experiments  indicate  a reduction  in  spe- 
cific total  resistance  of  0.029  for  the  greater  sidehull 
stagger  of  —85  m at  a speed  of  40  knots.  This  rep- 
resents an  impressive  saving  and  it  vindicates  the 
unusual  and  novel  nature  of  our  trimaran  concept. 


Figure  10:  Contours  of  free-surface  elevation 

mapped  onto  the  the  free  surface.  Camera  is  located 
between  the  center  and  starboard  hulls,  looking  aft 


This  reduction  in  total  resistance  would  be  dif- 
ferent at  prototype  scale  because  of  the  difference  in 
the  Reynolds  number. 

Secondly,  it  is  heartening  to  note  that  the  the- 
ory predicts  the  nature  of  the  resistance  reduction 
observed  in  the  experiments.  However,  at  the  de- 
sign speed  of  40  knots,  the  prediction  suggests  a 
reduction  of  only  0.021. 

Finally  we  turn  to  Figure  11(b),  which  is  a 
plot  of  the  ratio  of  the  two  resistances  ARt/Rt > 
rather  than  a difference  which  was  plotted  in  Fig- 
ure 11(a).  With  this  form  of  plot,  the  relative  over- 
all agreement  between  theoretical  predictions  and 
experiments  is  actually  more  encouraging.  From  a 
practical  viewpoint,  we  see  that  there  can  be  a rel- 
ative reduction  in  total  resistance  of  up  to  43%,  at 
a speed  of  38  knots. 


Figure  11:  Change  in  resistance:  Configuration  2 
and  Configuration  3 WRT  Configuration  1 

6 Conclusions 


6.1  Current  Work 

This  paper  contains  the  results  of  calm- water  resis- 
tance prediction  for  a high-speed  sealift  trimaran. 
The  two  numerical  codes  that  were  used  are  a 
thin-ship  theory  and  the  commercial  CFD  software 
FLUENT. 

It  has  been  demonstrated  that  the  thin-ship 
approach,  as  exemplified  by  the  results  plotted  in 
Figure  7,  Figure  8,  and  Figure  11,  provide  remark- 
ably reasonable  estimates  of  the  total  resistance,  if 
one  simply  adds  the  theoretical  wave  resistance  Rw, 
the  hydrostatic  drag  Rp,  and  the  ITTC  1957  fric- 
tional resistance  Rp.  It  should  be  borne  in  mind 
that  computational  times  are  typically  one  second 
per  speed,  using  a computational  grid  of  60  panels 


longitudinally  and  20  panels  vertically.  Thus,  the 
CPU  time  is  six  of  seven  orders  of  magnitude  less 
than  that  required  for  the  CFD  computations.  Of 
course,  the  CFD  results  are  considerably  more  use- 
ful, because  they  provide  data  for  the  complete  flow 
field  and  would  be  beneficial  to  the  naval  architect 
in  many  ways. 

Figure  7 has  been  used  to  demonstrate  the 
effectiveness  of  employing  a frictional  form  factor 
fp  = 1 + fc,  which  has  been  calculated  through  an 
estimate  of  the  increased  subhull  surface  velocities 
due  to  the  subhull-proximity  effects,  as  already  de- 
tailed. The  method  developed  so  far  results  in  a 
value  of  fp  ss  1.20.  That  is,  it  is  suggested  that  the 
frictional  resistance  in  such  cases  of  closely  spaced 
subhulls  is  increased  by  about  by  20%.  The  use  of 
such  a factor  in  Figure  8 would  bring  the  thin-ship 
predictions  into  nearly  perfect  alignment  with  the 
full  CFD  computations,  without  the  need  for  any 
significant  additional  CPU  time. 

The  CFD  results  contained  in  this  paper  were 
calculated  with  a rudimentary  grid  generation  ap- 
proach that  reduces  set-up  and  solution  time.  A 
goal  of  the  project  which  supported  this  work  was  to 
evaluate  the  codes  for  transfer  to  a shipyard  and/or 
design  office.  Often  in  the  application  of  CFD,  the 
grid  generation  is  a frustrating  hurdle  for  users  to 
overcome,  so  the  approach  here  is  very  simple  to  ef- 
fect, yet  it  does  not  seem  to  sacrifice  accuracy.  The 
overall  domain  is  divided  into  near-field  and  far-field 
regions  that  are  filled  with  polyhedral  and  hexahe- 
dral  cells  respectively. 

The  polyhedral  finite-volume  cell  that  was 
used  drastically  reduces  cell  skewness  and  overall 
cell  count.  The  reduced  skewness  causes  the  equa- 
tion coefficient  matrices  to  be  more  diagonally  dom- 
inant, and  thus  improves  solution  convergence  and 
robustness. 

The  accuracy  of  the  CFD  results  was  outstand- 
ing. The  error  relative  to  the  experimental  values 
was  on  average  less  than  2%  with  a maximum  of  4%. 
We  remind  the  reader  that  the  CFD  predictions  are 
done  using  the  experimental  values  of  sinkage  and 
trim. 

The  reconfigurable  catamaran  is  shown  to  dra- 
matically reduce  its  total  drag  by  extending  the  cen- 
terhull.  The  experiments  dictate  a 43%  reduction 
in  total  drag  at  model  scale.  This  reduction  is  suffi- 
cient to  encourage  further  investigation  of  the  con- 


cept,  namely  consideration  of  the  actuation  system 
for  the  extending  hull  and  requisite  gear,  and  the 
impact  on  the  hydrodynamic  performance.  Also, 
during  the  initial  design  period,  it  was  unclear  which 
sidehull  spacing  would  be  best  for  low  resistance. 
Interestingly,  it  was  the  narrowest  Configuration  6 
that  had  the  lowest  drag,  though  the  effect  was  very 
small. 


6.2  Future  Work 

The  authors  are  encouraged  by  the  quality  of  the 
results  from  both  the  thin-ship  theory  and  CFD. 
At  the  same  time,  there  are  issues  which  will  be 
valuable  to  investigate  further  to  expand  the  utility 
of  these  numerical  tools  for  the  designer. 

The  thin-ship  predictions  can  be  improved  by 
a more  careful  study  of  the  frictional  form  factor, 
which  is  thought  to  be  the  main  cause  of  the  total 
drag  estimate  being  a little  low.  The  reader  is  re- 
minded that,  in  the  literature,  there  is  a large  body 
of  evidence  indicating  that  for  some  vessels,  the  fric- 
tional form  factor  might  be  as  high  as  1.40  - far 
exceeding  the  model  value  of  around  1.20  indicated 
earlier. 

Another  point  is  the  estimate  of  the  length  of 
the  hollow,  which  is  assumed  to  grow  as  the  speed 
increases,  following  the  extensive  model  tests  re- 
ported by  Doctors  (2006a  and  2006b).  This  increase 
in  effective  hydrodynamic  length  of  the  vessel  low- 
ers the  predicted  wave  resistance.  Studies  of  this 
matter  are  also  in  the  planning  stage. 

The  CFD  calculations  are  very  expensive  in 
both  the  need  of  parallel  computational  hardware 
and  the  long  computational  time  to  arrive  at  the 
converged  solution.  Any  research  that  aims  at  ren- 
dering the  solution  algorithms  more  efficient,  or  re- 
ducing the  number  of  unknowns  would  be  valuable. 
It  is  natural  to  expect  the  solution  quality  to  im- 
prove with  increased  grid  density,  though  this  has 
not  been  proven  unequivocally  for  CFD  surface-ship 
resistance  predictions.  With  the  new  polyhedra  el- 
ement now  emerging  in  the  CFD  community,  and 
the  always  evolving  naval  combatant  hull  forms,  a 
detailed  grid-refinement  study  is  justified  as  a con- 
tinuation of  this  work,  with  the  optimism  that  less 
rather  than  more  dense  and  complex  grids  will  be 
required  to  arrive  at  useful  predictions. 
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